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We study Si:P donor electron spin decoherence due to anisotropic hyperfine (AHF) interaction 
with the surrounding nuclear spin bath. In particular, we clarify the electron spin echo envelope 
modulation (ESEEM) in the Si:P system and the resonance-like contributions from nuclear spins in 
various shells away from the P atoms. We suggest an approach to minimize AHF-induced decoher- 
ence by avoiding the resonances and orienting an applied magnetic field along directions that can 
periodically eliminate contributions from the dominant nearest neighbor atoms. Our remarkable 
agreement with experiment demonstrates nearly complete understanding of electron spin decoher- 
ence in Si:P when combining ESEEM, spectral diffusion, instantaneous diffusion, and spin-lattice 
relaxation. 
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I. INTRODUCTION 

The prospect of scalable solid state quantum infor- 
mation processing has created extensive interest in the 
study of coherent manipulation of single electron spins 
confined in semiconductor quantum dots or donor states. 
The clarification of single electron spin properties such 
as decoherence and quantum control should also be valu- 
able scientifically in the field of nanostructure physics and 
very useful technologically in the context of spintronics 
and spin quantum memory. 

Among the many semiconductor host materials, silicon 
is particularly enticing^ due to its deep-rooted connection 
to the modern microelectronics industry. Furthermore, Si 
provides a remarkably quiet environment for electron 2 ^ 
and nuclear— spin qubits since 28 Si, its most abundant 
isotope, has no net nuclear spin, and spin-orbit interac- 
tion is weak in Si. In this context, Si, with P donor elec- 
tron spins as qubits, is an ideal semiconductor material 
for quantum information processing. The donor electron 
spin dephasing time Ti as measured directly (and cal- 
culated theoretically) by spin echo decay is extremely 
long (many milliseconds in contrast to microseconds in 
GaAs) and can be further enhanced via isotopic purifi- 
cation. However, a strong anisotropic hyperfine (AHF) 
interaction in Si presents a formidable challenge for Si 
spin quantum computation; it produces finite spin de- 
phasing relatively quickly that cannot be eliminated by 
spin refocusing techniques, and could potentially nullify 
the advantage of long spin Ti time in Si. 

This paper details the effects of the AHF interaction 
on donor electron spin qubits in Si:P. We discuss this 
interaction and present our Hamiltonian formalism in 
Sec. HU In Sees. IIIII and IIV1 we formulate the problem 
of electron spin evolution and decoherence due to the 
AHF interaction in the context of free induction decay 
and (Halm) spin echoes, respectively. We compare the- 
oretical computations with experimentally observed spin 
echo envelopes in Sec. |V1 where we demonstrate nearly 



complete theoretical understanding of the decoherence 
of the Si:P electron spin qubit. Combining AHF-induced 
electron spin echo envelope modulation (ESEEM) with 
decoherence from spectral diffusion, instantaneous diffu- 
sion, and spin-lattice relaxation, we account for all major 
sources of decoherence and achieve nearly perfect quan- 
titative agreement with Hahn echo decay experiments. 
We believe that AHF interaction supplies the final piece 
of the puzzle with respect to understanding Si:P electron 
spin decoherence. Furthermore, Sec. I VII provides valu- 
able insight into AHF-induced decoherence and a pre- 
scription for its suppression. Without proper treatment, 
this decoherence can violate the stringent fault-tolerant 
requirements of qubit fidelity much sooner than the nom- 
inal T2 time suggests. While there are no existing spin 
echo techniques to generally remove AHF-induced echo 
modulations^ we suggest a concrete method for sup- 
pressing this decoherence and clarify requirements on the 
applied magnetic field for spin quantum computation in 
Si. We give concluding remarks in Sec. IVIII 



II. ANISOTROPIC HYPERFINE 
INTERACTION 

The hyperfine (HF) interaction between an electron 
and a nuclear spin describes the magnetic dipolar cou- 
pling between the two spin species^ With S denoting 
the spin operator of the electron and I that of the nu- 
cleus, the HF Hamiltonian is given by TLhf = I ■ A ■ S, 
where the hyperfine tensor A is given by 
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with the electron position measured relative to the site of 
the nucleus, \I/ the electron wavefunction, and 7s and 7/ 
the gyromagnetic ratios of the electron and the nucleus. 
The first term of Eq. (fT]) is the isotropic Fermi contact 
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HF interaction. The second term is anisotropic. Which 
part of the interaction is more important depends on the 
electron wave function. For example, the GaAs conduc- 
tion band minimum occurs at the T-point of the Bril- 
louin zone, where the electron Bloch function is mostly 
atomic s-type, so that HF interaction in GaAs between 
an electron near the conduction band minimum and the 
surrounding nuclear spins is essentially isotropic. On the 
other hand, the degenerate conduction band minimum 
for Si occurs close to the X-point of the Brillouin zone, 
where the electron Bloch functions have significant con- 
tributions from p- and d-atomic-orbitals^^ so that HF 
interaction between an electron near the conduction band 
minimum, such as an electron confined to a donor or 
a quantum dot, and the surrounding nuclear spins has 
strong anisotropic characteristics. Indeed, AHF interac- 
tion has been studied extensively in the Si:P system in 
the 1960s and 1970si±&±i The strength of AHF has been 
accurately measured^ and calculated^ for the phospho- 
rus donor electron. In the context of solid state spin 
quantum computation, however, much of the existing lit- 
erature only take into account the contact HF [first term 
in Eq. (fT])] in considering electron spin decoherence. 

To analyze how AHF leads to spin decoherence, we 
consider a single P donor in Si, with the donor-bound 
electron interacting with the P and randomly distributed 
29 Si nuclear spins. We assume the limit of a strong mag- 
netic field (> 100 mT is sufficient) applied in the z direc- 
tion such that electron spin flips are suppressed due to 
its large Zeeman energy. Since 7s ^> 7/, it is appropriate 
to take the limit where S z is conserved but not I z (of 
any nucleus). In this limit we write the Hamiltonian (in 
H = 1 unit) as H = Ho + ^2 n H n with 

Ho = us§*+Ap§,I?-w P I?, (2) 
Hn = A n S z I nz + B n S z I nx ' — uiil nz . (3) 

We separate the Hamiltonian into Ho , involving the elec- 
tron Zeeman energy and the donor nucleus, and H n , in- 
volving the nth 29 Si nucleus in the surrounding lattice 
(other Si isotopes have zero spin). In our notation, S, 
I p , and J n denote spin operators of the electron, P nu- 
cleus, and the nth 29 Si nucleus, respectively. I nx > gives 
the nuclear spin operator in an x'-axis orientation so that 
there is no S z I ny i contribution (thus we generally have 
a different x' orientation for each n). Given an applied 
magnetic field strength of B, we define lu\j = 7D-B as 
the Zeeman frequency for the electron, P nucleus, or a 
29 Si nucleus with □ = S,P, or I respectively. Ap de- 
notes the HF coupling between the electron and the P 
nucleus. Both contact HF and the S Z I Z part of the AHF 
interaction are contained in A n . The remaining AHF in- 
teraction in our strong field limit is contained in B n and 
gives the relevant anisotropy, mixing different directional 
components of S and J. 

Qualitatively, due to the anisotropic term B n S z I nx i in 
H n , the quantization axis for the precession of the 29 Si 
nuclear spin is dependent on the state of the electron 



spin. Conversely, the electron spin is affected by the pre- 
cession of the nuclear spin. The resulting electron spin 
free induction decay (FID) in Si:P has been explored in 



Ref. 12j and will be briefly reviewed in Sec. IIIII It is 



shown 12 that the donor electron spin could lose more 
than 1% of its coherence after only about 10 fis if a 29 Si 
atom is in one of the nearest neighbor (E-shell) sites. This 
would be disastrous for quantum computation where the 
error rate must stay below typical fault tolerance require- 
ments of 10~ 6 — 10~ 4 . Fortunately, as we will show in 
Sec. IVI1 AHF-induced decoherence may be drastically 
suppressed by applying precisely timed pulses and a suf- 
ficiently strong magnetic field along special directions of 
high symmetry. 



III. EFFECT OF ANISOTROPIC HYPERFINE 
ON FREE INDUCTION DECAY 

H n of Eq. §5§ for different n commute with each other 
so that the free evolution operator, in the interaction 
picture, becomes the product of evolution operators for 
each n: 
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We may write this evolution operator in the form 
Uo(t) = U+(t)P 1+ U -(t)P L 



Uo{t) = e~ mt = e- inat I I e- m »*. (4) 
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where Pj- = |T)(TI an d P\, = are up and down pro- 

jection operators for the electron spin In other words, 
U^it) denotes the free evolution of the nuclei given an 
electron spin that is either up or down. The free evolu- 
tion operator may take the form of Eq. ([5]) because the 
Hamiltonian only involves S z and no other electron spin 
operator. For the same reason, only dephasing occurs 
for the electron spin and thus only the off-diagonal ele- 
ments of the electron spin's density matrix may evolve. 
Assuming unpolarized (random) nuclei, after time t the 
off-diagonal element of the electron spin's density matrix 
becomes^ 2 - 
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where Au n = uj n+ - cJ n _, w n = {uj n + + w„_)/2, with 



(6) 




(7) 



and a 2 + b n = 1 (details in Ref. [l2j). iv n ± are preces- 
sion frequencies for nucleus n. Nucleus n will precess at 
a frequency of uj n + or uj n - given an up or down electron 
spin, respectively. In general, if the electron spin is in 
a superposition of up and down states, the nuclear spin 
dynamics will contain both these frequencies. As men- 
tioned before, numerical evaluation of Eq. © shows 12 
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FIG. 1: (Color online) (a) Maximum modulation depth 
[Eq. (|14[) ] in natural Si averaged over isotopic configurations 
with an applied magnetic field, B, parallel to the [001] lattice 
direction considering all shells (provided in Ref. 7[), j us ^ E- 
shell sites (nearest neighbors of the P donor), and all shells 
except the E-shell. This maximum depth gives the worst-case 
scenario of constructively interfering nuclear ESEEM con- 
tributions and is useful for identifying the major sources of 
modulations. When one shell of nuclei dominate, the worst- 
case maximum modulation depth is close to the actual ob- 
served modulation depth. Near the cancellation condition, 
<jJi ~ An/2, for each shell of nuclei is a peak labeled by the 
shell letter where that shell dominates; away from these peaks, 
the E-shell is seen to dominate, (b) Enlargement of the A- 
shell peak marking three field strengths, a, /3, and \, used 
by experiments presented in this Article, (c) Corresponding 
experimental Hahn echo decay at field strengths a and (3 for 
the same Si:P sample. Relatively high doping, 10 16 P/cm 3 , 
results in fast exponential relaxation due to instantaneous dif- 
fusion, but ESEEM is still observed. 



that the donor electron spin could lose more than 1% of 
its coherence in ~ 10 fis if a 29 Si atom is in one of the 
nearest neighbor (E-shell) sites. 



IV. ELECTRON SPIN ECHO ENVELOPE 
MODULATIONS 

The key question now is whether the AHF-induced 
electron spin decoherence can be suppressed. It is well 
known that spin echo techniques such as Hahn echo can 
be used to remove dephasing caused by the spatial vari- 
ation of local magnetic fields (the inhomogeneous broad- 
ening). However, the AHF-induced FID is a dynam- 
ical effect and, as such, cannot be removed by Hahn 
echo. Instead, AHF causes the echo envelope to oscillate, 
which is known within the electron spin resonance com- 
munity as ESEEM. 6 ' 13 i 14 This effect is particularly useful 
in chemistry for the identification of nuclear spin species 
of molecular sites&i^. The focus of the present study is 
to investigate ESEEM in the Si:P systemi&i&iZjiS and 
explore possible ways to reduce the decoherence effect of 
AHF interaction with 29 Si in the context of spin quantum 
computation. 



In the Hahn echo sequence, after initializing the elec- 
tron spin (which is necessary for experimental decoher- 
ence measurements, but is neither applicable to qubit 
preservation nor necessary in the current discussion), the 
system evolves freely for a time r, then a resonant pulse 
rotates the spin by 7r about an axis perpendicular to the 
magnetic field, then an echo is observed after another r 
in time (for shorthand, we use r — > n — > r to denote this 
sequence). For convenience in our formalism, though it 
makes no difference, we apply a second 7r pulse in order 
to bring the spin back to its original orientation (apart 
from decoherence): r — > 7r — > r — » n. Assuming an ideal 
applied n pulse and free evolution in the diagonal form 
of Eq. ([5]), the evolution over the echo sequence has the 
diagonal form 
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Mathematically, spin echo is the magnitude of the ex- 
pectation value of the electron spin (a measure of the 
fidelity of the qubit) after the quantum system evolves ac- 
cording to the echo sequence [Eq. ©]. As it is physically 
measured, the echo is the average spin of an ensemble 
of electron spins that undergo the Hahn echo evolution. 
Assuming only S z electron spin interactions, as we do, 
maximum decoherence occurs for an electron that is ini- 
tially perpendicular to the z direction (the direction of 
the applied magnetic field), and the spin echo is more 
precisely defined with such an initial electron spin state. 
In the experimental measurement, the ensemble of elec- 
tron spins are initialized to point in the same direction 
perpendicular to this axis. Where (...) is the quantum 
mechanical average over initial bath states, this Hahn 
echo envelope may be expressed as^ 9 - 



V(r) = 



u- 



(10) 



In computing the expectation value, we average over all 
initial nuclear states, assuming a completely disordered 
and unpolarized initial bath as in Sec. IIIII 

Electron spin echo decays as a function of r due to in- 
ternuclear interactions. In the current study we do not 
focus on this process, known as spectral diffusion, which 
has been studied by two of usi^ previously. On top of 
this decay, AHF coupling [Eq. ([3])] produces envelope 
modulations or ESEEMi&is Neglecting the internuclear 
interactions, the ESEEM due to each nucleus factors into 
the Hahn echo evolution separately [Eq. Q]. Given that 
only a fraction / of the Si nuclei have nonzero spin ( 29 Si), 
we have the following ESEEM amplitude averaged with 
respect to isotope configurations! 1 ^ 20 



V(r) = Y[[0-f)+fV n {r)] 
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1 - — (1 - cos (lo„+t)) (1 - cos (w„_t)) ,(12) 
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where tv n ± are nuclear precessional frequencies defined 
by Eq. 0. In the literature, k n is called the modula- 
tion depth parameter.— The maximum modulation (de- 
viation from I) of V n (r) is 2fc„, so that k n is a measure 
of modulation amplitude. In the "worst-case" scenario, 
when modulations from all nuclei combine constructively, 
the maximum possible modulation depth averaged over 
isotopic configurations is given by 

Max (1 - V(t)) = 1 [] t 1 2 / fc »] ■ ( 14 ) 

n 

In Fig. [1] (a), we show this maximum modulation as a 
function of field strength due to various nuclear shells 
(symmetry-related sets of lattice sitesi). We use exper- 
imentally determined contact and AHF coupling con- 
stants for 22 nuclear shells (which include about 150 
symmetry-related nuclear sites) taken from Ref. 0] and 
Ref. [IH. 

V. UNDERSTANDING EXPERIMENTAL 
OBSERVATION OF SPIN ECHO DECAY AND 
MODULATION 

The Hahn echo decay of Si:P donor electron spin has 
been studied in several recent experiments ! 15 i 16 ' 17 i 18 A 
previous quantum theory of nuclear-induced spectral dif- 
fusion (SD) by two of usi£ shows very good quantitative 
and qualitative agreement with these experimentally de- 
termined decay curves. However, the agreement is not 
perfect, especially at short times, because we previously 
neglected the effects of AHF interactions that produce 
modulations of the echo envelope, i.e., ESEEM. We now 
add this last piece to the puzzle to acheive truly re- 
markable agreement with experimental echo decay curves 
and thus demonstrate a nearly complete understanding 
of Si:P donor electron spin decoherence. 

Figure [21(c) shows excellent agreement of our ESEEM 
calculations with experimental data reported in Ref. 16]. 
The theory calculations do use five separate fitting pa- 
rameters: normalization, strain distribution width, re- 
laxation time, SD time, and a SD exponent. The first 
three of these parameters (described momentarily) may 
be fixed for all different directions of the applied mag- 
netic field. Thus, in Fig. [31 which shows comparison with 
experiment^ for ten different magnetic field directions, 
we use only two fitting parameters per curve. These 
two fitting parameters characterize the SD decay and are 
compared with the results of our SD theory^ in Fig. 2) 
Considering that we use, in our SD analysis, the Kohn- 
Luttinger envelope function within the effective mass ap- 
proximation for the donor electron,— which is known to 
be less reliable near the phosphorus atom , 22 ' 23 the agree- 
ment is quite good. There appears to be, however, some 
discrepancy between the fit and the theory for the SD 
exponent when the applied field is close along the [001] 
lattice direction: theory expects n = 2.3 and the fit yields 
n = 2.5. It is probably not coincidental that the nearest 



4 



B II [001] Magnetic Field Direction (9 = 0°) 




x (ms) 

FIG. 2: (Color online) AHF-induced ESEEM in Si:P with 
an applied magnetic field in the [001] direction, (a) Pure 
AHF-induced ESEEM for a single electron spin. The green 
"blob" is one curve with high frequency components, (b) 
Before matching the ESEEM to experiment, we must ac- 
count for strain effects in the ensemble of donor electrons 
(green); we must additionally sample at the same values of r 
as the experiment (red) yielding a stroboscopic effect, (c) 
Comparison with experiment (black). In our calculation, 
we combine decoherence effects of ESEEM, (non-Markovian) 
nuclear- induced SD, and (Markovian) exponential relaxation 
by simply multiplying them together. The orange curve 
gives ESEEM of our theory [red curve in (b)] multiplied by 
exp [-2t/Tr] exp [- (2r/Tsc) n ], where T R , Tsd, and n are 
fitting parameters for the relaxation time, SD time, and SD 
exponent, respectively. 



neighbor dipolar coupling vanishes when the applied field 
points along the [001] direction. Perhaps we have over- 
looked some interactions that become important when 
the dipolar interaction is weak. Overall, however, we do 
demonstrate good theoretical understanding of Si:P elec- 
tron spin decoherence in Figs. [HO andHJ 

We now describe the fitting parameters in more de- 
tail. The AHF-induced ESEEM [Eq. (II)] for a single 
electron spin, using the experimentally determined cou- 
pling constants for 22 nuclear shells from Ref. 0] is shown 
in Fig. 0(a). Random defects such as dislocations cause 
strain effects (such as electron population shifts between 
different Si valleys) that slightly alter the coupling con- 
stants for different donors. Strain effects result in nar- 
row distributions for the values of HF coupling constants 
and/or Zeeman frequencies and effectively dampen the 
ESEEM signal for an ensemble of spins. This is shown in 
Fig. [2h) where, in order to fit the experimental results, 
we assume a Gaussian distribution for HF frequencies 
with a 0.4% widths We believe strain is the culprit of 
this distribution of HF interaction strengths because of 
two reasons. First, experimentalists observe strain ef- 
fects orders of magnitude larger than this when applying 
a small stress on their sample, so it is quite plausible 
and likely that random defects are generating the strain, 
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FIG. 3: (Color online) AHF-induced ESEEM in Si:P for ten 
different curves corresponding to ten different magnetic field 
angles ranging from the [OOf] to the [ffO] directions. The 
plots are shifted in order to distinguish each angle. All fits 
use the same normalization, strain distribution width (0.4%), 
and relaxation time (Tr = 2.17±0.02 ms) parameters. There 
are two fitting parameters per curve: the SD time, Tsd, and 
the SD exponent, n. These fitting parameters are compared 
with our SD theory^ in Fig. H 



consistent with our excellent fit i ' 2 i 26 Second, other ne- 
glected interactions, such as electrons interacting with 
nuclei at different donors and dipolar interaction between 
nuclear spins, are much too small to account for this ef- 
fect. If we do not assume a distribution of HF coupling 
constants, the echo modulation would continue with a 
constant amplitude [see Eq. (fT!?)) ] until eventually nuclear 
spin relaxation kicks in. This is obviously not what is ob- 
served experimentally ] 15 i 16 i 17 i 18 

In addition to strain, Fig. HJb) shows the stroboscopic 
effect that emerges when we sample the same values of 
r as those reported in the experiment. The theoretical 
(orange) curve in Fig. [2jc) shows the ESEEM result of 
Fig. fflb) multiplied by exp [-2t/T r ] exp [- {2 T /T SD ) n } 
to account for independent effects of Markovian relax- 
ation and non-Mar kovian nuclear- induced SD. After we 
normalize the signal strength as an additional fit (the 
experiment only gives the Hahn echo decay on a rela- 
tive scale based on the strength of the observed signal) , 
we obtain excellent agreement with the experimental re- 
sults [black curve in Fig. [SJjc)]. Again, we use a total of 
five fitting parameters in Fig. [21 however, we use only 
two fitting parameters per curve in Fig. [3] and these two 
SD parameters are compared with theoretical results via 
Ref.[l| in Fig. H 

Although nuclear-induced SD and ESEEM both result 
from interactions with the nuclear spin bath, we have 
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FIG. 4: (Color online) Comparison between the SD fitting 
parameters (black triangles connected with dotted lines) of 
Fig.|3]and the theoretical predictions (red triangles connected 
with dashed lines) of our SD 19 theory. Right [left] triangles 
correspond to n [Tsd]; sizes approximate fitting uncertainty. 
The experimental n fit deviates from theory (2.30 ±0.05) only 
at small angles, where nearest neighbor dipolar flip-flop inter- 
actions approach zero. 



confirmed that it is appropriate to treat them indepen- 
dently. Nuclear-induced SD results from very small con- 
tributions from many thousands of pairs of flip-flopping 
nuclei.— The weak nuclear dipolar coupling that is re- 
sponsible for SD is too small to have any significant im- 
pact on the ESEEM-contributing nuclei that have a much 
stronger HF coupling to the electron. To be rigorous, we 
have made explicit cluster expansion calculations 1 ^ that 
include the AHF interaction. ESEEM emerges from one- 
cluster contributions, and the two-cluster contributions 
are negligibly changed with the introduction of AHF 
coupling. This calculation confirms the independence of 
these two decoherence channels. 

The Markovian relaxation is dominated by instan- 
taneous diffusion (ID) that results from interactions 
between the resonant electron spin donors. 6 Having 
a concentration of 8 x 10 14 donors / cm 3 in these 
experiments^ with half of those at resonance with the 
applied pulses (determined by whether the P donor nu- 
cleus has an up or down spin), yields an ID time of 
Tin ~ 3 ms; 6 - accounting for most of the Tr ~ 2.2 ms 
relaxation. A further relaxation process with a decay 
time of about 8.3 ms would be needed to account for 
the remaining contribution to Tr (the inverse of the de- 
cay times are additive). With these experiments per- 
formed at 8 Kfi 6 - this may be attributed to temperature- 
dependent spin-lattice relaxation; 8.3 ms is fairly con- 
sistent with reported 10 — 25 ms temperature-dependent 
relaxation time in Si:P at 8 K (a wide uncertainty range 
is due to strong temperature dependence and offsets in 
temperature calibration) ^ 

It is important to emphasize here that the theoretical 
fitting for the echo modulations (separate from spectral 
diffusion and relaxation) uses experimentally obtained 
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AHF interaction strength. The only variable fitting pa- 
rameter for echo modulations is the 0.4% distribution of 
the AHF interaction strength, which allows us to repro- 
duce the decay of the echo modulation. In the meantime, 
as indicated in Fig. [4] the fit for spectral diffusion is an 
attempt to use a relatively simple functional form to rep- 
resent our theoretical calculations in Ref . . Therefore 
the seemingly large number of fitting parameters involved 
here is not a blind attempt to obtain the best possible fit 
to the experimental measurements. Instead, they mostly 
provide a more transparent representation of a more so- 
phisticated theory, and the remarkable agreement we ob- 
tain here illustrates that all the main decoherence mech- 
anisms are represented, most likely, in our description. 



VI. SUPPRESSING 
ANISOTROPIC-HYPERFINE-INDUCED ESEEM 

After building confidence in our theoretical approach 
by comparing our results with experiments, we now ad- 
dress the question of how to suppress spin decoherence 
induced by AHF interaction. More specifically, we will 
show that this AHF-induced decoherence may be dras- 
tically reduced by carefully choosing the strength and 
direction of the applied magnetic field and by applying 
spin echo pulses with precisely prescribed timing. 

One interesting and important feature of the maxi- 
mum modulations shown in Fig. [T] (a) is that a peak 
occurs when u>i ~ A n /2 (with A n positive) for each shell 
of atoms. At each such cancellation condition, as it is 
dubbed, the Zeeman and HF energies of nucleus n can- 
cel when the electron spin is up but not down, freeing 
the nuclear spin from conservation of energy constraints 
conditional upon the state of the electron spin. Mathe- 
matically, uj n+ is minimized [Eq. I)?])] so that k n is at (or 
very near) its maximum, resulting in modulation depth 
peaks. This effect is shown experimentally by comparing 
the two echo decay curves in Fig. [T] (c); the curve cor- 
responding to a field strength closer to the center of the 
A-shell peak clearly exhibits stronger echo modulations. 
It turns out that the experiments discussed in SeafVlwere 
performed with a magnetic field strength denoted as x m 
Fig. Q] (b) , which is very close to the center of the A-shell 
peak, where echo modulations are particularly strong. 

It is clear from this discussion that to minimize AHF- 
induced decoherence, cancellation conditions for all the 
shells with finite AHF coupling constant should be care- 
fully avoided by properly selecting the applied magnetic 
field strength (or, in electron spin resonance, the cor- 
responding microwave cavity frequency). Furthermore, 
away from the cancellation condition peaks, the E-shell 
nuclei (nearest neighbors to the P nucleus) have the 
strongest AHF coupling by far, so that they dominate 
the echo modulations by more than an order of magni- 
tude, as seen in Fig. [T] (a). Remarkably, the echo mod- 
ulation due to these dominating E-shell nuclei can be 
effectively removed at special magnetic field orientations 




FIG. 5: (Color online) Special applied magnetic field direc- 
tions that allow effective removal of echo modulation contri- 
butions due to E-shell nuclei (the four nearest neighbors to 
the P donor). The arrows and translucent sheets, respectively, 
indicate directions from the P atom parallel and perpendicu- 
lar to the applied magnetic field. Sites in these direction give 
no anisotropic contribution (B n = 0) because these directions 
are, by symmetry, along the principal axes of the A tensor of 
Eq. Q; thus, the top site in (b) and the two in-plane sites in 
(c) do not contribute to ESEEM. In each of the three cases, 
the E-shell sites that do contribute are magnetically equiva- 
lent. 



(Fig. [5| , if we exploit the periodic restoration of electron 
spin coherence in the presence of nuclei that are seen by 
the electron as magnetically equivalent. This restoration 
arises because V n (r) = 1 [Eq. I|12p] when r is a multi- 
ple of 2ir/w n ± (either + or — ). Note that such periodic 
restoration does not generally occur in the free induction 
case of Sec. IHII 

For special magnetic field orientations shown in Fig. 
the contributing E-shell sites are magnetically equivalent 
with the same {w n _, uj n +}', thus, the electron spin is pe- 
riodically restored at the same values of r regardless of 
isotopic ( 29 Si) configuration. In this way, E-shell con- 
tributions can effectively be eliminated as exemplified in 
Fig. [6l By orienting the magnetic field in one of the var- 
ious special directions, the effects of all E-shell nuclei are 
simultaneously eliminated at periodic values of r. 

To understand the periodic restoration of ESEEM in 
the presence of magnetically equivalent nuclei, note that 
Uo(t) simply generates precession for each nucleus at a 
frequency of to n ±, respectively. For instance, given an 
electron spin that is up, the nuclear spin returns to its 
initial orientation after waiting for a time that is a multi- 
ple of 2it/lu + . That is, assuming all nuclei being consid- 
ered are magnetically equivalent, U^{2irm/ uj±) = 1 for 
any integer m. If we consider a Hahn echo sequence with 
t = 2tt/w + , we have [see Eq. ©] U± aha = U^(t)U±(t) = 
Uq(t). Thus the evolution of the magnetically equiv- 
alent nuclei is independent of the electron spin, since 
UH ilhn = ^Hah„j so that the electron qubit is fully de- 
coupled from these nuclei (apart from the effects of the 
neglected nuclei that are not magnetically equivalent to 
the others). This property is common to any balanced 
sequence in which an initially up or down electron (or the 
separate components of a superposition state) spends an 
equal amount of time being up and down. Specifically, we 
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FIG. 6: (Color online) Echo modulations, 1 — V(t), in nat- 
ural Si with an applied field of 1 T in the [001] direction 
corresponding to Fig. [5] (a). When sampling r at multiples 
of 2n/to+ (or 2n/u>-), the E-shell nuclei give no contribution 
to the echo modulations. Error bars, vertically asymmetric 
because of the logarithmic scale, correspond to the standard 
deviation resulting from random isotopic configurations; those 
with down arrows extend below the visible range. 



10 kHZ. This implies that we are limited to t C 100 fis. 
On the other hand, the situation in a different experiment 
or sample could be quite different and would depend on 
the precise effects of strain on the E-shell nuclei which is 
outside of the scope of this paper. Also, depending on the 
architecture of the proposed quantum computer, it may 
be possible to calibrate pulse timing to individual qubits 
rather than an ensemble, circumventing the uncertainty 
caused by strain. 

Errors in the applied magnetic field direction cannot 
be fixed with calibration. If the magnetic field direction 
deviates from one of the three special directions in Fig. [SI 
the contributing E-shell nuclei will no longer be equiva- 
lent and nuclei that are supposed to be inert [e.g., the 
top nucleus in Fig. [5] (b) and the two leftmost nuclei in 
Fig. [5j (c)] will, in general, contribute to ESEEM. In or- 
der to address the question of sensitivity to the magnetic 
field direction, we first consider how angular uncertainty 
affects the uncertainty of io n ± : 



1 



^n± 



- { ±A n /2-u Jl ) — + -B n — 



may apply our timing trick to eliminate E-shell modula- 
tions for the Carr-Purcell sequence o 27 ' 28 or concatenated 
dynamical decoupling sequence a 29 i 30 that can be much 
more effective than the Hahn echo for mitigating other 
types of decoherence such as spectral diffusionj 28 i 30 

It is important to determine how sensitive our tech- 
nique of suppressing AHF-induccd decoherence is to er- 
rors in timing and magnetic field direction. Because the 
E-shell contributions oscillate at frequencies of oj n - and 
uj n+ , the timing of the pulse sequence in order to remove 
this effect must be accurate with errors that are small 
compared with the period of these oscillations, 27r/w„_ 
and 2ir/oj n+ . For magnetic fields that are > IT, these 
oscillation frequencies [Eq. Q] are dominated by the nu- 
clear Zeeman frequency, uj n ± rs luj, so that the periods 
are approximately 2ir/u>i ~ 100 ns, as is apparent in 
Fig. E 

The timing is also affected by systematic uncertainty 
in either co n+ or (whichever we choose for synchro- 
nization) in a way that worsens in time and number of 
oscillation periods. The resulting uncertainty in time is 
given by 



5t = 



( SuJni 



(15) 



Then, the requirement that 5t <C 2-K/uj n ± yields the fol- 
lowing restriction: 



r < 



2tt 



(16) 



When dealing with an ensemble, the uncertainty of to n ± 
will be subject to strain effects. In Sec. [Vj we empiri- 
cally determined strain effects causing a 0.4% distribu- 
tion, which equates to a frequency uncertainty of about 



We again assume that u n ± ~ loj, so that 



max 



5uj n± 


i 


SA n 


se 


w — max 


59 


2 n 



(17) 



This is valid when max„|A„| > max n |B n |, so that the 
5B n contribution is negligible. We have verified this nu- 
merically by calculating how A„ and B n vary with re- 
spect to changing the magnetic field direction about an 
axis that rotates from [100] to [011]. Although we did 
not analyze full two-dimensional changes in direction, we 
do extract a characteristic scale of 50 kHz/deg for both 
5A n /S9 and SB n /S9 of the E-shell nuclei, which we ex- 
pect to be general. In addition to w n ±, we must also con- 
sider how an error in the magnetic field direction away 
from [111] or [011] affects k n for those nuclei which are 
not expected, ideally, to contribute. For such nuclei, 



5 2 k n 



k r , 



B n =0 



I5 2 k n 



2 59 2 
2u 2 



(59)\ 



(18) 



B„=0 



'S BJ59 



SB r , 



(19) 



We find that k n B/(59) 2 ~ 0.5% T 2 / deg 2 . Thus, with 
a one degree error in the magnetic field direction and 
a field strength of 1 T, the sites that are not supposed 
to contribute to ESEEM when the field is oriented along 
[100] or [011] could, if occupied by a 29 Si nucleus, actually 
provide modulations up to about 1%. This, however, 
scales quadratically with the angle error, so an error of 
1% of a degree would keep these errors below the 10~ 6 
threshold often quoted for quantum error correction. 
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VII. CONCLUSION 

We have studied Si:P donor electron spin decoherence 
due to AHF interaction, which is an important dephasing 
mechanism in Si. We clarify the electron spin echo enve- 
lope modulation in the Si:P system and the resonance- 
like contributions from nuclear spins in various shells 
away from the P atoms. Our theory is in excellent quanti- 
tative agreement with experiment. Most importantly, we 
suggest an approach to minimize the decoherence effect of 
AHF interaction by avoiding the cancellation conditions 
and orienting an applied magnetic field along directions 
that can periodically eliminate the contributions from the 
dominant E-shcll nuclei. We quantitatively analyze the 



sensitivity of this strategy to errors in timing and mag- 
netic field orientation. 

In principle, one can eliminate the problem of nuclear 
spin decoherence through isotopic purification, remov- 
ing the 29 Si. After purification, those donors that still 
have 29 Si nuclei in the E-shell may be disqualified as 
qubits. Our strategy provides an alternative solution to 
this problem that provides the quantum computer archi- 
tect with more flexibility. 

We wish to thank David Cory, Steve Lyon, Sergei 
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especially would like to thank Alexei Tyryshkin for his 
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